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DESCRIPTION 



Motion artifact compensation 

The present invention relates to the field of image processing, for example,, in medical 
5 applications. In particular, the present invention relates to a method of motion artifact 
compensation in a projection data set, to data processing devices and to respective 
computer programs. 

In a CT scanner with a polychromatic source of radiation, such as a polychromatic x-ray 
1 0 source, a polychromatic x-ray beam passes through matter and photons are absorbed or 
scattered according to the properties of the matter they are passing through. The 
introduction of cone-beam computer tomography (CT) systems offers a couple of 
benefits. Compared with a single slice CT system, the time for data acquisition can be 
reduced, the output of the x-ray tube is used more efficiently, leading to a simpler heat 
15 management of the tube, and there is no longer a need for compromising z-resolution in 
order to achieve an acceptable scan time. However, a major problem using cone-beam 
CT systems is to find a proper reconstruction algorithm. A large number of publications 
within the last few years has addressed the problem of reconstruction of cone-beam 
projections acquired during a helical path of the source-detector system. Thus far, exact 
20 or quasi-exact algorithms exist only for a PI and a 3-PI-acquisition. Alternatively, an 
approximate algorithm maybe used, such as PI-FBP, WEDGE or WEDGE-PI. 

Motion is one of the most critical sources of artifacts in helical cone-beam CT. This is 
especially true if state of the art exact methods are used for reconstruction. Typically, 
25 these methods use an angular range of exactly n-PI of projection data for the 

reconstruction of each object point. This feature implies that the first and last ray used 
for reconstruction of a certain object point are exactly the same, but taken from opposite 
sides. Motion artifacts appear as arcs, which are due to an inconsistency between the 
first and the last ray. Approximate methods, like WEDGE mitigate this problem by so- 
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called over-scan weighting. In this case, a trade-off between good suppression of motion 
artifacts by using a large over-scan range and best image quality for the static case by 
using a small over-scan range has to be found. 

It is an object of the present invention to provide for an improved motion artifact 
compensation. 

In accordance with an exemplary embodiment of the present invention as set forth in 
claim 1, the above object may be solved by a method of motion artifact compensation in 
a projection data set of an object of interest, wherein the projection data set is acquired 
by means of a source of electromagnetic radiation generating a beam and by means of a 
radiation detector detecting the beam. According to this exemplary embodiment of the 
present invention, a difference between a first ray and a second ray is determined, after 
which the projection data set is compensated for a motion artifact on the basis of the 
difference, resulting in a motion artifact compensated projection data set. After that, the 
object of interest is reconstructed from the motion artifact compensated projection data 
set, resulting in a motion artifact compensated image. The first ray and the second ray 
correspond to projection data of the projection data set. The first ray and the second ray 
are opposite rays and the difference between the first ray and the second ray is due to the 
motion artifacts. 

In other words, an image is taken from an object of interest and projection data, which 
has been measured twice (for the first time by the first ray and for the second time by 
the second opposing ray) are compared to each other. If a difference is determined 
between the data resulting from the first ray and the second ray, the projection data set is 
compensated for a motion artifact, depending on the difference. After motion artifact 
compensation, a reconstruction of the object of interest is performed on the basis of the 
motion artifact compensated projection data set. 
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Advantageously, comparing the projection data resulting from the first ray and the 
projection data resulting from the second ray may be performed automatically during 
the scan, which may lead to a fast motion artifact compensation. 

5 According to another exemplary embodiment of the present invention as set forth in 
claim 2, the determination of a difference between the first ray and the second ray 
further comprises the steps of selecting the first ray and the second ray on the basis of 
the projection data and determining whether the difference between the first ray and the 
second ray is greater than a predetermined threshold, wherein, if the difference is greater 
10 than the predetermined threshold, a motion artifact compensation of the projection data 
set is performed. 

Advantageously, according to this exemplary embodiment of the present invention, no 
motion artifact compensation of the projection data set is performed if the difference 
15 between the first ray and the second ray does not exceed the predetermined threshold 
value. Therefore, by selecting a certain threshold value, the sensitivity of the motion 
artifact compensation algorithm maybe set according to individual needs. 

Another exemplary embodiment of the present invention is set forth in claim 3, wherein 
20 the second ray is interpolated from adjacent rays. Therefore, even if there is no original 
second ray which is exactly opposite to the first ray, such a second ray may be generated 
by interpolation from adjacent rays, which may be rays corresponding to projection data 
acquired in the neighborhood of the interpolated projection data corresponding to the 
interpolated second ray. 

25 

According to another exemplary embodiment of the present invention as set forth in 
claim 4, the object of interest comprises a plurality of object points, wherein a 
reconstruction of a first object point of the plurality of object points is performed by an 
exact reconstruction algorithm. Furthermore, if the motion artifact results from a motion 
30 of the first object point, the motion artifact is compensated for by a low pass filtering of 
the projection data in the region of the motion artifact before the reconstruction of the 
first object point by the exact reconstruction algorithm. 
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Advantageously, according to this exemplary embodiment of the present invention, this 
may smear out inconsistency over a large area and may thus result in a less severe 
motion artifact. 

5 

According to another exemplary embodiment of the present invention as set forth in 
claim 5, the source of radiation moves around the object of interest and the exact 
reconstruction algorithm uses projection data resulting from one of half a revolution and 
three half revolutions of the source of radiation. 

10 

Advantageously, by using projection data resulting from half a revolution of the source 
of radiation, opposing first and second rays may be easily and quickly detected, even 
with small detector arrays. 

15 Another exemplary embodiment of the present invention is set forth in claim 6, wherein 
characteristics of the low-pass filtering correspond to properties of the projection data in 
the region of the motion artifact. Advantageously, according to this exemplary 
embodiment of the present invention, a filter threshold may be dynamically selected, 
depending on the severity of the motion artifact. 

20 

According to another exemplary embodiment of the present invention as set forth in 
claim 7, the object of interest comprises a plurality of object points, wherein a 
reconstruction of a first object point is performed by an approximate reconstruction 
algorithm. An over-scan range is used for reconstruction of the first object point. If the 
25 motion artifact results from a motion of the first object point, the motion artifact is 
compensated for by increasing the over-scan range. Advantageously, this procedure 
allows for a local optimization of the image quality. 

According to another exemplary embodiment of the present invention as set forth in 
30 claim 8, the first object point belongs to a Pi-line, on which motion has been detected 
and the increase of the over-scan range corresponds to properties of the projection data 
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in the region of the motion artifact. Advantageously, according to this exemplary 
embodiment of the present invention, the over-scan range maybe set according to the 
severity of the motion artifact. In other words, a rather dominant motion artifact may 
lead to a large over-scan range and a smaller motion artifact may lead to a small over- 
5 scan range. This dynamic setting of the over-scan range on the basis of, for example, the 
difference between the first ray and the second ray, may lead to a fast and individual 
motion artifact compensation. 

According to another exemplary embodiment of the present invention as set forth in 
10 claim 9, the approximate reconstruction algorithm is one of a WEDGE algorithm and a 
Pi-filtered back projection algorithm (PI-FBP). Thus, a high sensitivity for motion 
artifacts may be achieved. 

According to another exemplary embodiment of the present invention as set forth in 
15 claim 10, the source of electromagnetic radiation is a polychromatic x-ray source, which 
moves along a helical path around the object of interest. Furthermore, the beam has one 
of a cone beam geometry and a fan beam geometry. 



The application of polychromatic x-rays is advantageous, since polychromatic x-rays are 
20 easy to generate and provide a good image resolution. Furthermore, since the geometry 
of the CT scanner system may be of different designs, such as, for example, cone beam 
or fan beam geometry, a method for an exemplary embodiment of the present invention 
may be applied to a plurality of different scanner systems and may not be limited to CT 
scanner systems. 

25 

According to another exemplary embodiment of the present invention as set forth in 
claim 1 1, a data processing device is provided, which comprises a memory for storing a 
data set and a data processor for performing motion artifact compensation in a 
projection data set of an object of interest, wherein the data processor is adapted for 
30 performing the following operation: loading the data set acquired by means of a rotating 
source of electromagnetic radiation generating a beam and by means of a radiation 
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detector detecting the beam; determining a difference between a first ray and a second 
ray; compensating the projection data set for a motion artifact on the basis of the 
difference, resulting in a motion artifact compensated projection data set; and 
reconstructing the object of interest from the motion artifact compensated projection 
5 data set, resulting in a motion artifact compensated image. The first ray and the second 
ray correspond to projection data of the projection data set and are opposite rays. 
Furthermore, the difference between the first ray and the second ray is due to the motion 
artifact. 

10 Advantageously, this may allow for improved image quality of images comprising 
motion artifacts. 

The present invention also relates to a computer program, which may, for example, be 
executed on a processor, such as an image processor. Such computer programs may be 

15 part of, for example, a CT scanner system. The computer programs, according to an 
exemplary embodiment of the present invention, are set forth in claim 12. These 
computer programs may be preferably loaded into working memories of data 
processors. The data processors are thus equipped to carry out exemplary embodiments 
of the methods of the present invention. The computer programs may be stored on a 

20 computer readable medium, such as a CD-ROM. The computer programs may also be 
presented over a network such as the WorldWideWeb, and may be downloaded into the 
working memory of a data processor from such networks. 

It may be seen as the gist of an exemplary embodiment of the present invention that 
25 motion artifacts in a projection data set are detected by determining a difference 
between a first ray and a second ray, which are opposite to each other and therefore 
correspond to projection data, which do not show a difference in the absence of motion. 
After determining the difference, the projection data set is compensated for a motion 
artifact depending on the magnitude of the difference. This may provide for a better 
30 image quality in reconstructed images. 
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These and other aspects of the present invention will become apparent from and 
elucidated with reference to the embodiments described hereinafter. 

Exemplary embodiments of the present invention will be described in the following, 
5 with reference to the following drawings: 

Fig. 1 shows a simplified schematic representation of an embodiment of a computed 
tomography (CT) scanner according to the present invention. 

10 Fig. 2a shows a schematic representation of a virtual planar detector according to a cone 
beam reconstruction. 

Fig. 2b shows a parallel beam geometry according to a PI-FBP algorithm. 
15 Fig. 3 shows a geometry for back-projection according to a PI-FBP algorithm. 
Fig. 4 shows a geometry for determining s according to a PI-FBB algorithm. 
Fig. 5 shows a schematic representation of a virtual focus-centered detector. 

20 

Fig. 6 shows the geometry for determining V w according to a WEDGE algorithm. 

Fig. 7 shows an example of a weighting function w (©). 

25 Fig. 8 shows a histogram of a Pl-line-deviation for a moving (dynamic) and non-moving 
(static) heart. 

Fig. 9 shows Pi-line deviations for a moving heart (left image) and a non-moving heart 
(right image) in the form of an axial slice. 

30 

Fig. 10 a shows a reconstructed axial slice for a moving heart, which is not compensated 
for motion artifacts according to the present invention. 
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Fig. 10 b shows two reconstructed axial slices with the PI + plus + EPSELON algorithm 
for the moving heart. 

Fig. 1 1 shows a schematic representation of a Pi-line. 

5 

Fig. 12 shows a schematic representation of a 3 Pi-line. 

Fig. 13 shows a flow-chart of an exemplary embodiment of method of motion artifact 
compensation according to the present invention, 

10 

Fig. 14 shows a flow-chart of another exemplary embodiment of a method of motion 
artifact compensation according to the present invention. 

Fig. 15 shows an exemplary embodiment of an image processing device according to 
15 the present invention for executing an exemplary embodiment of a method in 
accordance with the present invention. 

Fig. 1 shows an exemplary embodiment of a CT (computed tomography) scanner 
system according to the present invention. With reference to this exemplary 
20 embodiment, the present invention will be described for the application in medical 
imaging. However, it should be noted that the present invention is not limited to the 
application in the field of medical imaging, but maybe used in applications such as 
baggage inspection to detect hazardous materials, such as explosives, in items of 
baggage or other industrial applications, such as material testing. 

25 

The scanner depicted in Fig. 1 is a cone beam CT scanner. The CT scanner depicted in 
Fig. 1 comprises a gantry 1, which is rotatable around a rotational axis 2. The gantry is 
driven by means of a motor 3. Reference numeral 4 designates a source of radiation 
such as an x-ray source, which, according to an aspect of the present invention, emits a 
30 polychromatic radiation. 



-9- 



PHDE040044 EPP 



Reference numeral 5 designates an aperture system which forms the radiation beam 
emitted from the radiation source to a cone shaped radiation beam 6. 

The cone beam 6 is directed such that it penetrates an object of interest 7 arranged in the 
5 center of the gantry 1, i.e. in an examination region of the CT scanner and impinges 
onto the detector 8. As may be taken from Fig. 1, the detector 8 is arranged on the 
gantry 1 opposite the source of radiation 4, such that the surface of the detector 8 is 
covered by the cone beam 6. The detector 8 depicted in Fig. 1 comprises a plurality of 
detector elements. 

10 

During a scan of the object of interest 7, the source of radiation 4, the aperture system 5 
and detector 8 are rotated along gantry 1 in the direction indicated by arrow 16. For 
rotation of the gantry 1 with the source of radiation 4, the aperture system 5 and the 
detector 8, the motor 3 is connected to a motor control unit 17, which is connected to a 
1 5 calculation unit 1 8. 

In Fig. 1, the object of interest is disposed on a conveyor belt 19. During the scan of the 
object of interest 7, while the gantry 1 rotates around the patient 7, the conveyor belt 19 
displaces the object of interest 7 along a direction parallel to the rotational axis 2 of the 

20 gantry 1 . By this, the object of interest 7 is scanned along a helical scan path. The 

conveyor belt 19 may also be stopped during the scans to thereby measure single slices. 
Instead of providing a conveyor belt 19, for example, in medical applications, where the 
object of interest 7 is a patient, a movable table is used However, it should be noted 
that in all of the described cases it is also possible to perform a circular scan, where 

25 there is no displacement in a direction parallel to the rotational axis 2, but only the 
rotation of the gantry 1 around the rotational axis 2. 

The detector 8 is connected to the calculation unit 18. The calculation unit 18 receives 
the detection result, i.e. the read-outs from the detector element of the detector 8, and 
30 determines a scanning result on the basis of the read-outs. The detector elements of the 
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detector 8 may be adapted to measure the attenuation caused to the cone beam 6 by the 
object of interest. Furthermore, the calculation unit 18 communicates with the motor 
control unit 17 in order to coordinate the movement of the gantry 1 with motor 3 and 20 
or with the conveyor belt 19. 

5 

The calculation unit 1 8 may be adapted for reconstructing an image from read-outs of 
the detector 8. The image generated by the calculation unit 1 8 may be output to a 
display (not shown in Fig. 1) via an interface 22. 

10 The calculation unit, which may be realized by a data processor, may also be adapted to 
perform a motion artifact correction compensation in the image based on the read-outs 
from the detector elements of the detector 8. According to an aspect of the present 
invention, this motion artifact compensation may be performed by determining a 
difference between a first ray and a second ray, wherein the first ray and the second ray 

15 are opposite rays and wherein the difference between the first ray and the second ray is 
due to a motion artifact. In the following, the motion artifact may be compensated for 
on the basis of the determined difference, resulting in a motion artifact compensated 
projection data set. After that, the object of interest may be reconstructed from the 
motion artifact compensated projection data set, resulting in a motion artifact 

20 compensated image. 

Furthermore, the calculation units may be adapted for performing a motion artifact 
compensation in the projection data sets of the object of interest by performing the 
following operation: loading the data set acquired by means of a rotating source of 

25 electromagnetic radiation generating a beam and by means of a radiation detector 
detecting the beam. After that, a difference between a first ray and a second ray is 
determined, wherein first and second rays are opposite rays and wherein the difference 
between the first and the second ray is due to the motion artifact. First ray and second 
ray correspond to projection data of the projection data set. After that, the projection 

30 data set is compensated for a motion artifact on the basis of the difference, resulting in a 
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motion artifact compensated projection data set. Following the compensation, a 
reconstruction of the object of interest from the motion artifact compensated projection 
data set may be performed, resulting in a motion artifact compensated image, 

5 Furthermore, as may be taken from Fig. 1, the calculation unit 1 8 may be connected to a 
loudspeaker 21 to, for example, automatically output an alarm. 

In the following, a short description of the basic formalism of several approximate 
helical cone-beam reconstruction algorithms is presented. The X-ray focus a(X) is 
1 0 moving along a helical trajectory: 

a(X) — {-R -cos/l, -RsinA,, h/l + z 0 ) 

The Pitch Pis then 

15 

P = 27ih. 

One possible detector is a virtual planar detector which is always perpendicular to the 
central ray and moves with the focus in the z-direction. It can be characterized for each 
20 focus-position X with two coordinates, u and v which are defined as follows (see Fig. 
2a): 

u = R- tan y v = z-h<Z-z 0 

25 This means that each projection value is characterized by the triple (u, v, A,). The 
virtual detector is spanned by two normal vectors e u and e v (pointing in the positive z- 
direction and is therefore e z ). The vector e f maybe defined according to 
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10 



which is perpendicular to the virtual detector and points to the X-ray focus. Another 
virtual detector is the Pi-window (or Tam-Danielsson window) which is delimited by 
two successive turns of the helix. It can be expressed by the coordinates of the virtual 
planar detector according to: 



^ Lower — V < V upper 3 



with 



2 

Vu P per,P.-window = h • (1 + ^) • - arctan^-) 
v Lower, pi-window = h • (1 + ^5-) • (— — - arctan — ) 

The Pi-window guarantees a data acquisition scheme, which is in principle complete for 
15 an exact reconstruction of the object function, but also non-redundant. 

The original PI- window may be generalized to the so called n-PI-window by a 
symmetric extension of the PI- window by one or more times the helical pitch. The 
upper and lower boundaries of the n-PI-window expressed by the coordinates of the 
20 virtual planar detector are: 

v Up per,n-Pi-window = h - (1 + - (— - arctan-) 

V Lower, n-PI-Window - « " U + "j^f) * V ^~ ~ arCtan— ) 

25 The importance of the PI and the n-PI window is due to the fact that from all projection 
data measured at a certain source position X, only the data, which are located on the 
lower or upper boundary of these windows are also measured from the opposite 
direction. Thus, only these data are used in the present invention to detect motion. 
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The PI-FBP Algorithm 

The PI-Filtered back-projection Algorithm (PI-FBP Algorithm) is an approximate 
5 reconstruction algorithm for a helical trajectory. It works in principle like that: 

• it operates only on projection data in the Pi-window 

• the data are rebinned into parallel beam geometry 

• pre- weighting of the rebinned data with the cosine of the cone-angle (this 

10 compensates the different path lengths of the X-rays and guarantees that objects 
which are homogeneous in z-direction can be exactly reconstructed) 

• ramp filtering row by row 

• back-projection of the filtered projection data 

15 The rebinning of the data corresponds to a change in the coordinate systems, expressed 



by 



(u, v, X)<^(t, e, s) 



20 with (see Fig. 2b) 




—u 



0 = X 4- y = A - arctan(— ) 



s:=z-h-0-z o = 



v 





25 
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10 



15 



In the above formulas, y is the fan angle, t is the distance of the line integral to the z- 
axis when projected into the x-y-plane and 9 is the angle between this line-integral and 
the X-axis. A short calculation shows that the s-coordinates of the upper and lower 
boundaries of the n-PI-window are: 



S Upper,n-Pl-wIndow ~ S Lower, n-PI- window ~ 1Z~ (Equation (1)) 



with h;—k] for the n-PI-window. 

2 2 

As was already mentioned, after rebinning into parallel beam geometry, the projection 
data are weighted with the cosine of the cone angle k which is: 



cos k = 



Afterwards, the data are ramp-filtered in the t-direction (as in 2D-parallel beam 
reconstruction). At the end, for each voxel with the coordinates (x, y, z) a back- 
projection is done to compute the object function g(x,y,z): 

g(x, = — P Filtered, Parallel V> % y, Z, )) 

Q n=n Sunrise {x,y,z) 



The following relations are given (see Fig. 3): 



y = arctan( ) 

R +x • cos X + y ' siaA 

y = arcsm( — ) 
R 

Q = X + y 
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t = — x • sin 0 + y • cos 6 

The distance p, from (x,y,z) to the point Q is: 
5 //(x ? 7 ? #) = cos <9 + y - sin # 

For the s-coordinate, the following relations are given (see Fig. 4): 

Az = ^R 2 -t 2 ±/j(x,y,0) 
tenet ^R 2 ~t 2 

10 

Az = z — hX — ^ 0 , 
s = z-h0-z Q =Az Det -hy, 

15 and therefore: 

z {z-hX-zA^R 2 -t 2 
s = -hy+ . v ^ . 

^R 2 -t 2 +x-cos0 + y-sin& 

For every voxel, the back-projection may only be performed for x-rays belonging to the 
20 Pi-window which corresponds to an angular range of AG = n . This means that one has 
to calculate for every voxel when the X-ray (so called Pi-line) enters the Pi-window (for 

this Pi-line we have s = -^~) and determine its coordinates (t, 0, s) SR where SR 

stands for sunrise. The same Pi-line exists for the "sunset" (X-ray leaves the PI- 
window) , with 9 SS = G SR + k, t ss = -t SR and s ss = -h%l 2 . 



25 
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The WEDGE-Algorithm 

The WEDGE-Algorithm is also an approximate reconstruction algorithm for a helical 
trajectory. It works in principle like that: 

5 

• the data are rebinned into parallel beam geometry (but different when compared 
with the PI-FBP algorithm, see below) 

• pre-weighting of the rebinned data with the cosine of the cone-angle (this 
compensates the different path lengths of the X-rays and guarantees that objects 

10 which are homogeneous in z-direction can be exactly reconstructed) 

• ramp filtering row by row 

• back-projection of the filtered projection data 

The rebinning of the data into parallel beam geometry is identical with the PI-FBP 
1 5 algorithm. However, it is different with respect to the coordinate describing the detector 
height. For the PI-FBP algorithm we already have shown the coordinate transformation 
from a virtual planar detector to the Pi-Geometry: 

(u 9 v, A) <=>(*, 0, s) 

20 

For the WEDGE-algorithm, the corresponding coordinate transformation is: 

(u, v, A)<=>(t, 0, v w ). 

25 The remaining transformation v w (u, v, X) may be derived from the following relations 
(see Fig. 5 and Fig. 6): 



v w =Az = z — hA, — z 0 
Az' = z' - hX - z 0 
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20 



Az' Vr^?" 

v w =(z'-/*0-z o ) * 



The following relation between s in the PI-FBP geometry and v w in the WEDGE 
geometry is obtained: 

R 

5 S J R2 _ f2 ~ v » (Equation (2)) 

As for the PI-FBP algorithm, after rebinning into the WEDGE beam geometry, there is 
a pre-weighting of the projection data with the cosine of the cone angle to guarantee an 
exact reconstruction for objects, which are homogeneous in the z-direction. Afterwards, 
the ramp filtering in t-direction is applied. At the end, the back-projection is done for 

10 every voxel (x,y,z) with the filtered projection data. For the WEDGE algorithm, the 
back-projection is done over the whole detector height, not just the PI- window. This 
means that the angular range AO is usually larger than % . As an immediate 
consequence, there must be a weighting function w(0) which is used during the back- 
projection process and which compensates for the fact that some rays are twice present 

1 5 and others only once (see Fig. 7). However, one must be aware that such a weighting 
function maybe problematic since two rays which differ in 0 by % are not really 
parallel to each other. As a consequence, for objects with strong gradients in z-direction, 
the WEDGE algorithm usually shows more artifacts in the reconstructed image than the 
PI-FBP algorithm. 



WEDGE-PI Algorithm 



The WEDGE-PI algorithm is nearly identical to the WEDGE- Algorithm, but uses only 
projection data in the Pi-window. The difference with respect to the PI-FBP method is 
25 the filtering. 
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As already mentioned, the projection data in the Pi-window are in principle sufficient 
for an exact reconstruction of the measured object, but also non-redundant. This means 
that in principle no weighting of the projection data has to be done which compensates 
for the fact that some line integrals are measured more than once. 

5 

However, reconstruction algorithms which operate solely on the PI- window are much 
more sensitive to motion artifacts. During the back-projection process, the first and the 
last ray (Pi-line of the "sunrise" and Pi-line of the "sunset") are identical and should 
contain the same projection values. This is not true anymore, when motion is 
10 happening. This inconsistency in the projection values immediately manifests itself in 
streaking artifacts in the reconstructed image. 

Algorithms which perform the back-projection over an angular range larger than n are 
much less sensitive to motion artifacts. The reason is that the inconsistencies in the 
15 projection data set are somehow smeared out onto a larger areas in the reconstructed 
image and therefore not as easily visible (no distinct streaks). 

In order to improve image quality it may be advisable to use a "larger" back-projection 
window (larger than the PI- window) in image areas where motion has happened, but to 
20 use a relatively "small" angular range for the back-projection (Pi-window or at least not 
substantially larger than Pi-window) when the motion is "not so strong". A motion- 
detecting algorithm and subsequent adaptation of the angular range for the back- 
projection (WEDGE-PI+Epsilon Algorithm) according to an exemplary embodiment of 
the present invention comprises the following steps: 

25 

• for every voxel (x,y,z) the Pl-lines for sunrise and sunset are determined 

• if the projection values for the Pl-lines are (nearly) identical, one takes the PI- 
window for back-projection 

• if the projection values for the Pl-lines are rather different, one takes the full angular 
30 range (over the full detector height) for back-projection. 
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In the following, the WEDGE-PREpsilon Algorithm will be described in more detail. 
The algorithm works on projection data in WEDGE-geometry. 

Ideally, for a non-moving object the projection values of the Pi-Lines of every voxel 
5 should be the same: 

p(t,0,v w (s = nh/2)) = p(r-t,0 + ?r,v w (s = -nhlTj) 

(see equations (1) and (2)). In reality, they are not the same and one can define the 
relative deviation of the projection values of the Pi-Lines of every voxel: 

d 9 m= \p(t, 0, v w (s = 7rh/ 2)) - pirU 0 + 7T y V w (s=-7Zh/ 2))| 
10 Re/ * mzx(p(t, 0,v w (s = 7ih/ 2)), pi-t, 0 + tt 9 v w (s = -xh 1 2) ■ 

In the following, this quantity is called "Pi-Line-Deviation". 

The program, according to an exemplary embodiment of the present invention, may 
15 now loop over all voxels which have to be reconstructed. As already mentioned above, 
for each voxel the values of t, 0, and v w corresponding to the sunrise of the Pi-window 
are calculated. Calculating a projection value at these coordinates may generally require 
a three-dimensional interpolation since the real data are sampled discretely in all 
coordinates. In the following embodiment, a tri-linear interpolation is used. 

20 

The angular range for the back-projection of every voxel is adjusted individually by 
choosing two thresholds x 1 and x 2 : 

• if d Rel (x, y, z) < x 1 , one chooses A0(x, y, z) = % + s 1 
25 • if d Re ,(x, y, z) > x 2 , one chooses AG(x, y, z) = % + s 2 

• if t 1 < d Rel (x,y,z) < x 2 , one chooses 
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A0(x, y,z) = n + s i + 7* el X \ } • (s 2 - s,) . 

As an alternative, % + s 2 may be replaced by the angular range provided by the WEDGE 
algorithm (the default "Overscan" window). 

5 

For the thorax phantom a histogram depicting the relative Pi-Line deviations for the 
moving and non-moving heart is shown in Fig. 8. The motion-induced Pi-Line 
deviations are between 0.3% and 2%. Therefore, as already mentioned above, a careful 
determination of the Pi-Lines and their corresponding projection values may be 
10 mandatory so that the "signal" (d Re , ) may be separated from "discretization noise". Fig. 
9 shows the PI-Line-deviations for one axial slice for the moving heart (left image) and 
non-moving heart (right image). The effect of the moving heart is obvious. 

Fig. 10 a shows a reconstructed axial slice for a moving heart, which is not compensated 
1 5 for motion artifacts according to the present invention. On the left side of the picture, e.g., 
stripes resulting from a motion of parts of the heart can clearly be identified. Furthermore, a 
pronounced shadow artifact near the vertebra is present. 



Fig. 10b shows two reconstructed axial slices for a moving heart, compensated for motion 
20 artifacts according to an exemplary embodiment of a method according to the present 
invention, where s 1 was chosen as 30° and the angular range was adjusted in the interval 

[n + 8 1 ; Overscan] with the parameters x 1 =0.0025 and x 2 = 0.02 . The motion artifacts 
are clearly suppressed (left image of Fig. 10 b); furthermore, there is no pronounced 
shadow artifact near the vertebra present (right image of Fig. 10 b). 

25 

Fig. 11 shows a schematic PI measurement geometry of a CT scanner system for 
performing the method according to an exemplary embodiment of the present invention. 
For an exact PI reconstruction algorithm, a detector area is used, which corresponds to 
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the detector area depicted in Fig. 11, fitting between one cycle of the radiation source, 
which moves along a helical path around the object of interest and which emits a cone 
beam of polychromatic radiation. Pi-lines are located at the upper and lower boundaries 
of the detector 30. At position 34, the source 4 emits a beam, which comprises a first 

5 ray 32, striking a certain object point 33 and hitting detector 30 at position 35. After half 
a revolution of the source 4 along the helical path 3 1 to position 35, the source 4 emits a 
cone beam, which comprises a second ray 36, running from position 35 of source 4, past 
object points 33 to position 34, which is now located on the detector 30. Therefore, the 
rays 32 and 36, which are emitted at different times from different sides of the helical 

10 path are true opposite arrays. 

Fig. 12 shows a schematic 3 -PI measurement geometry for a CT scanner system for 
performing a method according to an exemplary embodiment of the present invention. 
The radiation source 4 moves along helical path 3 1 and emits a cone beam of 

15 electromagnetic radiation. At position 34, the source 4 emits a cone beam which 
comprises the first ray 32, passing the object point 33 and hitting the detector 30 at 
position 35. After one and a half revolutions of the source 4 along the helical path 31, 
the source is positioned at point 35 and emits a cone beam, comprising the second ray 
36, which passes the object point 33 and hits the detector 30 at position 34. Since the 

20 detector 30 moves along the helical path 31 and is always opposite radiation source 4, 
the detector 30 is located on the opposite side of its position depicted in Fig. 12, when 
source 4 has been moved to position 35. Therefore, the first ray 32 and the second ray 
36 are opposite rays, passing the same object point 33 and may therefore be used for 
motion detection and compensation of a motion artifact arising from a motion of object 

25 point 33 according to a method of an exemplary embodiment of the present invention. 

Fig. 13 shows a flow-chart of an exemplary embodiment of a method of motion artifact 
compensation according to the present invention. After the start in step SI, an 
acquisition of the projection data set is performed in step S2. After that, the acquired 
30 data set is loaded into the memory of a data processor and a first ray and a second ray 
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are selected in step S3, wherein the first ray and the second ray are selected on the basis 
of the projection data. After that, in step S4, it is determined whether the difference 
between the first ray and the second ray is greater than a predetermined threshold value. 
The threshold value may be set by a user or it may be set from the software side, 
5 depending on properties of the projection data. If the difference between the first ray 
and the second ray is not significant, meaning that it is less than or equal to the 
predetermined threshold value, no motion artifact correction is performed (step S5). On 
the other hand, if the difference between the first ray and the second ray is greater than 
the predetermined threshold value, a motion artifact compensation is performed in step 

10 S6 by low-pass filtering. Advantageously, the first ray or the second ray may be 

interpolated from adjacent rays, for example, if no original second ray may be found 
which is opposite to the first ray. In other words, if a motion artifact results from a 
motion of an object point which is passed by the first ray and the second ray, the motion 
artifact may be compensated for by a low-pass filtering of the projection data in the 

1 5 region of the motion artifact before the reconstruction of the first object point by an 
exact reconstruction algorithm. Then, in step S7, a reconstruction of the object point is 
performed by the exact reconstruction algorithm, which may use, according to an 
exemplary embodiment of the present invention, projection data resulting from half a 
revolution (PI) or three half revolutions (3 -PI) of the radiation source. 

20 

Furthermore, the characteristics of the low-pass filtering, such as filter strength and 
frequency dependence, may, according to an exemplary embodiment of the present 
invention, correspond to properties of the projection data in the region of the motion 
artifact. In other words, the low-pass filtering may be performed in regions where 
25 motion has been detected in an adaptive manner. 

The motion artifact compensation method ends in step S8. 

Fig. 14 shows a flow-chart of another exemplary embodiment of a method of motion 
30 artifact compensation according to the present invention. The method of motion artifact 
compensation starts with step SI. The acquisition of a projection data set of an object of 
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interest by means of a source of electromagnetic polychromatic radiation generating a 
cone beam and by means of a radiation detector detecting the cone beam is performed in 
step S2. After that, in step S3, a first ray and a second ray are selected, wherein the first 
ray and the second ray correspond to projection data of the projection data set and 
5 wherein the first ray and the second ray are opposite rays. After that, in step S4, the first 
ray and the second ray are compared to each other and it is determined whether there is 
a significant difference between the first ray and the second ray. If there is no significant 
difference between the two opposing rays, no motion artifact correction is performed 
(step S5). On the other hand, if a significant difference between the first ray and the 
10 second ray has been detected, the over-scan range is increased, leading to a motion 

artifact compensation in the data corresponding to the first ray and the second ray. This 
is performed in step S6. 

It should be noted that the object of interest, which is scanned by the CT scanner 
15 system, comprises a plurality of object points. One of these object points, the first object 
point, is passed by both the first ray and the second ray. The reconstruction of the first 
object point is performed by an approximate reconstruction algorithm, wherein an over- 
scan range is used for reconstruction of the first object point. By increasing the over- 
scan range of the object point, a suppression of motion artifacts may be increased. 
20 Therefore, if a difference between the first ray and the second ray exceeds a certain 
predetermined threshold value, the corresponding projection data is considered to 
comprise motion artifacts resulting from a motion of the first object point. In 
consequence, the motion artifact is compensated for by increasing the over-scan range 
of the first object point. Advantageously, the amount of the increase of the over-scan 
25 range corresponds to properties of the projection data in the region of the motion 

artifact. Therefore, the over-scan range used for reconstruction of an object point may 
be increased adaptively, depending on the properties of the projection data. This 
procedure may allow for a local optimization of the image quality. 



30 



-24- 



PHDE040044 EPP 



After that, in step S7, a reconstruction of the object of interest is performed by an 
approximate reconstruction algorithm, such as a WEDGE algorithm, a WEDGE-PI 
algorithm, a WEDGE-PI+Epsilon algorithm or a Pl-filtered-back projection algorithm, 
which are described above in more detail, 

5 

The method ends in step S8. 

Fig. 15 depicts an exemplary embodiment of an image processing device according to 
the present invention, for executing an exemplary embodiment of a method in 

10 accordance with the present invention. The image processing device depicted in Fig. 15 
comprises a central processing unit (CPU) or image processor 151 connected to a 
memory 152 for storing an image depicting an object of interest, such as a patient. The 
image processor 151 may be connected to a plurality of input/output network or 
diagnosis devices, such as an MR device or a CT device. The image processor is 

15 furthermore connected to a display device 154, for example, a computer monitor, for 
displaying information or image computed or adapted in the image processor 151. An 
operator may interact with the image processor 151 via a keyboard 155 and/or other 
output devices, which are not depicted in Fig. 15. 

20 Furthermore, via the bus system 153, it is also possible to connect the image processing 
and control processor 151 to, for example, a motion monitor, which monitors a motion 
of the object of interest. In case, for example, a lung of a patient is imaged, the motion 
sensor may be an exhalation sensor. In case the heart is imaged, the motion sensor may 
be an electrocardiogram (EGG). 

25 



-25- 



PHDE040044 EPP 



CLAIMS 



1. A method of motion artifact compensation in a projection data set of an object of 
interest, wherein the projection data set is acquired by means of a source of 
electromagnetic radiation generating a beam and by means of a radiation detector 
detecting the beam, the method comprising the steps of: 

determining a difference between a first ray and a second ray; 
compensating the projection data set for a motion artifact on the basis of the 
difference, resulting in a motion artifact compensated projection data set; 

- reconstructing the object of interest from the motion artifact compensated 
projection data set, resulting in a motion artifact compensated image; 

- wherein the first ray and the second ray correspond to projection data of the 
projection data set; 

wherein the first ray and the second ray are opposite rays; and 

- wherein the difference between the first ray and the second ray is due to the motion 
artifact. 

2. The method according to claim 1, 

- wherein the determination of a difference between the first ray and the second ray 
further comprises the steps of: 

selecting the first ray and the second ray on the basis of the projection data; 
determining, whether the difference between the first ray and the second ray is 
bigger than a predetermined threshold; 

- wherein, if the difference is bigger than the predetermined thresholds a motion 
artifact compensation of the projection data set is performed. 
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3. The method according to claim 1, wherein the second ray is interpolated from 
adjacent rays. 

4. The method according to claim 1, 

- wherein the object of interest comprises a plurality of object points; 

- wherein a reconstruction of a first object point of the plurality of object points is 
performed by an exact reconstruction algorithm; and 

- wherein, if the motion artifact results from a motion of the first object point, the 
motion artifact is compensated for by a low pass filtering of the projection data in 
the region of the motion artifact before the reconstruction of the first object point 
by the exact reconstruction algorithm. 

5. The method according to claim 4, 

- wherein the source of radiation moves around the object of interest; and 

- wherein the exact reconstruction algorithm uses projection data resulting from one 
of half a revolution and three half revolutions of the source of radiation. 

6. The method according to claim 4, wherein characteristics of the low pass filtering 
correspond to properties of the projection data in the region of the motion artifact. 

7. The method according to claim 1, 

- wherein the object of interest comprises a plurality of object points; 

- wherein a reconstruction of a first object point is performed by an approximate 
reconstruction algorithm; 

- wherein an over-scan range is used for reconstruction of the first object point; and 

- wherein, if the motion artifact results from a motion of the first object point, the 
motion artifact is compensated for by increasing the over-scan range. 
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8. The method according to claim 7, 

- wherein the first object point belongs to a Pi-line on which motion has been 
detected; and 

- wherein the increase of the over-scan range corresponds to properties of the 
5 projection data in the region of the motion artifact. 

9. The method according to claim 7, wherein the approximate reconstruction algorithm 
is one of a WEDGE algorithm and a Pi-filtered back-projection algorithm. 



10 10. The method according to claim 1, 

- wherein the source of electromagnetic radiation is a polychromatic x-ray source; 

- wherein the source moves along a helical path around the object of interest; and 
wherein the beam has one of a cone beam and a fan beam geometry. 

15 1 1 . A data processing device comprising: 

- a memory for storing a data set; 

- a data processor for performing motion artifact compensation in a projection data 
set of an object of interest, wherein the data processor is adapted for performing the 
following operation: 

20 - loading the data set acquired by means of a rotating source of electromagnetic 
radiation generating a beam and by means of a radiation detector detecting the 
beam; 

- determining a difference between a first ray and a second ray; 

- compensating the projection data set for a motion artifact on the basis of the 
25 difference, resulting in a motion artifact compensated projection data set; 

- reconstructing the object of interest from the motion artifact compensated 
projection data set, resulting in a motion artifact compensated image; 

- wherein the first ray and the second ray correspond to projection data of the 
projection data set; 
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- wherein the first ray and the second ray are opposite rays; and 

- wherein the difference between the first ray and the second ray is due to the motion 
artifact. 

5 12. A computer program for performing motion artifact compensation in a projection 
data set of an object of interest, wherein the computer program causes a processor to 
perform the following operation when the computer program is executed on the 
processor: 

- loading the data set acquired by means of a rotating source of electromagnetic 
1 0 radiation generating a beam and by means of a radiation detector detecting the 

beam; 

determining a difference between a first ray and a second ray; 

- compensating the projection data set for a motion artifact on the basis of the 
difference, resulting in a motion artifact compensated projection data set; 

15 - reconstructing the object of interest from the motion artifact compensated 
projection data set, resulting in a motion artifact compensated image; 

- wherein the first ray and the second ray correspond to projection data of the 
projection data set; 

- wherein the first ray and the second ray are opposite rays; and 

20 - wherein the difference between the first ray and the second ray is due to the motion 
artifact. 
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ABSTRACT 

Motion artifact compensation 

Motion is one of the most critical sources of artifacts in helical cone-beam CT. By 
comparing opposite rays corresponding to projection data, the amount of motion may be 
estimated and, in the following suppression of corresponding motion artifacts may be 
performed according to an exemplary embodiment of the present invention. The method 
of motion artifact compensation may be implemented in both approximate 
reconstruction algorithms and exact reconstruction algorithms. Advantageously, motion 
during the data acquisition is detected automatically and related motion artifacts may be 
suppressed adaptively. 
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